Transient photon production in a QGP 
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' Abstract. We discuss the shortcomings of a formula that has been used in the Hterature to compute 

I the number of photons emitted by a hot or dense system during a finite time, and show that the 

transient effects it predicts for the photon rate are unphysical. 
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Tij- , Electromagnetic radiation has long been thought to be a good probe of the early 

^ I stages of heavy ion collisions In recent years, a new "real-time" approach has been 

proposed in order to compute out-of-equilibrium effects for photon production in a dense 
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O . equilibrated quark-gluon system, which originate in its finite life-time [ 2] . The result was 

^ I that there are important transient effects that make the yield much larger than what would 

have been expected by simply multiplying the equilibrium rates by the corresponding 



X 



amount of time. This unexpectedly large photon yield was the starting point of many 
discussions regarding the validity of this approach [3^, JJ. 

In this note we show in a systematic and simple way that the expression of the 
^ . photon yield obtained in [2] relies on unphysical assumptions. We show that the standard 

canonical formalism in the 5-matrix approach leads to this expression, provided the 
electromagnetic interactions are unduly turned on and off at finite initial and final times. 
The simplicity of this derivation allows us to exhibit the illegitimate character of the 
expression used to predict the transient effects (For details, see Ref. 

Let us consider a system of quarks and gluons, and denote by Lagrangian. 
We couple the quarks to the electromagnetic field in order to study photon emission 
by this system, and denote ^e.m. the Lagrangian of the electromagnetic field, and 
the term that couples the quarks to the photons. The complete Lagrangian is therefore 
if = ^^^^ + ^e.m. +^<?r' Where ^^^^ = -(l/4)G;i^Gr +?(iVx -^^W -m) V^W ; 
^e.m. = -(l/4)i7iv^^^ ; -^qy = -e\f/{x)4^{x)\if{x) . Here and G^^ are respectively 
the gluon field and field strength, and F^v the photon field and field strength, and i// 
the quark field (only one flavor is considered here), g is the strong coupling constant, and 
e is the quark electric charge. We denote collectively by ^int the sum of all the interac- 
tion terms. The number of photons measured in the system at some late time is given by 
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the following formula 2(o{dN/d^xd\) = (l/y)Epoi. a (l/Z)Tr[p(/',-)a2^(k)«!,m (k)]. 
Here, V is the volume of the system and Z = Tr(p(?,)) the partition function. The 
sum runs over the physical polarization states of the photon. p{ti) is the density op- 
erator that defines the initial statistical ensemble. The "in" states and operators of the 
interaction picture are free, and are defined to coincide with those of the Heisenberg 
picture at the initial time t = tj. The "out" states and operators are those used in or- 
der to perform the measurement. In principle, the measurement should take place af- 
ter the photons have stopped interacting, i.e. one should count the photons at a time 
t +00 so that they are asymptotically free photons. Here, for the sake of the argu- 
ment, we are going to define the "out" states and fields at some finite time tf. This 
means that we assume that electromagnetic interactions have been turned off before 
the time tf, for this measurement to be meaningful. These "out" states and fields are 
related to the "in" states and fields by means of the "5-matrix": |aout) = 5'^|ctin) ; 

'^outC^) = ^^'^^m^)S ' ^iut^C^) = "^^^in , with an S matrix given in terms of 

the interaction as S = [/ {tf,ti) = jll d'^x^\ni{^m{x)), where denotes the path- 

ordering. Expressing the "out" creation and annihilation operators in terms of their "in" 
counterparts, and these in terms of the corresponding "in" fields, we arrive at: 

^""J^^^ ^ 4^)*(k)£W(k)(<9,,-f/a))(<9,„-/a)) 
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In this formula, the time path ^ goes from ti to tf along the real axis, and then back to 
ti, the time derivatives should be taken before xq and jq are set equal to tf, and (+) / (— ) 
stand for the upper/lower branch of ^. 

Now, one can perform an expansion in the electromagnetic coupling constant, while 
keeping strong interactions to all orders. This is motivated by the very different mag- 
nitude of the electromagnetic and the strong coupling constants. We can also sum over 
the photon polarizations, and combine the order and the order results. Then, for a 
system which does not contain any real photons initially, i.e. for which ny((u) = 0, the 
number of photons produced per unit time and per unit phase-space at the time tf is 

-\-oo 

dN 2 f dE sm{{E-(o){tf-ti)) ^^ 
^""did^ = J 2^ n^4£,k) . (2) 
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We have assumed that the initial density matrix describing the distribution of quarks 
and gluons is such that the photon polarization tensor, n^^(Mo, vo;k), is invariant under 
time translation. This formula is equivalent to the formula obtained by Boyanovsky et al 
10]. However, we have derived it here within the framework of an 5-matrix formulation. 
This was not the standard 5-matrix approach though, as some extra hypothesis and some 
extensions have been used. A first consequence of eq. Q is that it gives back the usual 
formula for the photon production rate at equilibrium if one takes to infinity the time 
tf at which the measurement is performed (which amounts to turn off adiabatically the 



electromagnetic interactions only at asymptotic times). At this point, the main question 
is whether the finite tf generalization of this formula makes sense as a photon production 

rate, as invoked in JS] . One may be tempted to interpret the number operator aout^out as 
the number operator at the time tf for photons still interacting with the system, but there 
is no warranty that this definition of the number of photons agrees with the number of 
photons as measured in a detector, precisely because they are not asymptotically free 
states. The only possibility to argue safely that this operator indeed counts observable 
photons is to assume that the system does not undergo electromagnetic interactions after 
the time tf. In any case, it is clearly unphysical to keep tf finite: either we are trying 
to measure non asymptotically free photons, or we have to turn off the interactions at 
a finite time tf. A similar problem arises at the initial time. The derivation we have 
used for eq. assumed that there is no dependence on e in the initial density operator 
p{ti), which is possible only if there are no electromagnetic interactions in the initial 
state. Moreover, imposing = in the initial state is also forbidding electromagnetic 
interactions before It is also known that eq. © is plagued by very serious pathologies 
that appear as ultraviolet divergences. Firstly, the r.h.s. of eq. ^ turns out to be infinite at 
any fixed photon energy ft) due to some unphysical vacuum contributions, i.e. processes 
where a photon is produced without any particle in the initial state. Secondly, the 
remaining terms, even if they give a finite photon production rate, lead to an energy 
dependence of this rate which is too hard for being integrable: one would conclude 
based on this formula that the total energy radiated as photons per unit time by a finite 
volume is infinite, which clearly violates energy conservation. It was claimed in fs*] 
that the vacuum terms could be discarded simply by subtracting to the r.h.s. of eq. Q 
the same formula evaluated in the vacuum. This indeed has the desired effect, but is 
a totally ad hoc prescription because nowhere in the derivation of the formula appears 
this subtraction term. Similarly, the authors of fs*] suggested that the divergence that 
appears in the total radiated energy can be subtracted by multiplying the creation and 
annihilation operators used in the definition of the number of photons by some wave 
function renormalization constants. However, no such constants appear in the derivation: 
the operators aoutj'^out be related to their "in" counterparts directly by means of the 
5- matrix. We conclude that the transient effects predicted in [2] are unphysical. 



ACKNOWLEDGMENTS 

The authors acknowledge CAPES/COFECUB for financial support. E.S.F. is partially 
supported by CNPq, FAPERJ and FUJB/UFRJ. 



REFERENCES 

1. F. Gelis, Nucl. Phys. A 715, 329 (2003). 

2. S.Y. Wang, D. Boyanovsky, Phys. Rev. D 63, 051702 (2001); S.Y. Wang, D. Boyanovsky, 
K.W. Ng, Nucl. Phys. A 699, 819 (2002). 

3. G.D. Moore, Communication at the CERN workshop on Hard probes in heavy ion colli- 
sions at LHC, March 11-15, 2002, Geneva, Switzerland. F. Arleo, et al, hep-ph/031 1 131^ 



J. Serreau, JHEP 0405, 078 (2004). 

E. Fraga, E Gelis and D. Schiff, hep-ph/03 12222 

D. Boyanovsky, H.J. de Vega, Phys. Rev. D 68, 065018 (2OO3); |Iiep^0M1156| 



